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ABSTRACT 


The  purpose  of  the  study  reported  upon  herein  is  to  investigate  the 
potential  existence  of  a  clustering  of  system  poles  and  teros  in  the  i -domain 
This  potential  is  shown  to  exist.  A  means  of  separating  the  clustering  is 
proposed  by  mapping  from  the  complex  z— domain  onto  the  cmnpleac  ^domain.  The 
foiKjn  for  so  doing  is  to  ameliorate  computational  difficulties  that  can 
arise  as  the  result  of  closely  clustered  roots. 
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SECTION  I.  INTRODUCTION. 

A.  Study  Objective. 

The  purpose  of  this  study  is  to  investigate  the  potential  existence  of 
a  clustering  of  system  poles  and  zeroe  in  the  z-d amain,  if  the  potential 
exists  *  a  mathematical  means  of  obviating  the  problem  is  to  be  determined* 

The  study  effort  is  composed  of  a  single  task  described  in  the  Scope  of 
Work  contained  in  Technical  Requirement  No*  T-0208  entitled  "Application  of 
Digital  Technology  to  Guidance  and  Control  Theory"  dated  21  May  1979**  The 
requirements  inposed  by  this  task  are  contained  in  the  referenced  document 
and  re-stated  (for  the  convenience  of  the  reader)  below. 

Task  I  requirements: 

"2.1  Investigate  to  determine  if  a  sensitivity  problem 
due  to  clustering  of  the  system  poles  and  zeros  exists  for 
the  parameter  space  technique  of  digital  system  analysis/design. 

2.2  If  a  sensitivity  problem  (i.e. ,  clustering  of  the  system 
poles  and  zeros)  is  found  in  the  investigation  of  Para.  2.1,  a 
technique  shall  be  developed,  to  mathematically  spread  apart  the 
system  poles  and  zeros  via  intermediate  complex  transform 
techniques  to  eliminate  the  sensitivity." 

B.  Study  Schedule. 

The  requirements  of  Task  I  were  completed  during  the  period 1  August 
1979  (contract  initiation)  through  30  September  1979.  They  are  reported  upon 
in  this  report. 

c.  tacUati  of.  Kmart. 

This  technical  report  covers  the  efforts  performed  in  completing  the 
requirements  of  Task  I  of  the  study.  As  stated  in  the  referenced  contract 
requirements:  "The  required  documentation  shall  be  a  technical  report." 
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SECTION  n.  PROBLEM  INVESTIGATION . 

As  Indicated  in  Section  2.1  of  the  referenced  Technical  Requirement ,  an 
Investigation  was  made  to  determine  if  a  sensitivity  problem  due  to  clustering 
of  the  system  poles  and  zeros  exists  far  the  parameter  space  technique  of  digital 
system  analysis/design.  As  will  be  shown  in  the  sequel*  any  technique  that 
utilizes  design  in  the  complex  z-domain  is  subject  to  this  problem.  If  this 

problem  is  shown  to  exist,  then  the  parameter  plane  technique  (among  others)  is 

2 

subject  to  it.  Hence  it  is  prudent  to  describe  what  the  nature  of  the  problem 
is  (Section  A,  below)  and  see  how  likely  it  is  to  occur  in  practice  (Section  B, 
below). 

a.  itotars  Pretelsn* 

The  nature  of  the  problem  can  be  seen  in  the  complex  z-domain.  Consider 
that  the  dynamics  of  a  system  are  represented  by  three  roots,  such  as  system  poles: 
two  that  are  complex  conjugates  and  one  that  is  real.  Further  suppose  that  they 
are  all  located  arbitrarily  close  to  the  point  z  -  +1.  This  situation  is  shown 
in  Fig.  1,  where  the  three  roots  are  designated  z^  ,  z2  *  and  -  *2  (the  over- 
bar  is  used  to  indicate  the  complex  conjugate  of  a  quantity): 

zt  “  d,  (2-1) 

and 

*2  -  pe1*  -  x  +  i  y,  (2-2) 

Where  d,  x,  and  y  are  real  quantities  and 

-  (  x?  ♦  /  )*,  (2-3) 

0  -  tan  "1  (y/x).  (2-4) 

The  complex  variable,  z,  is  defined  in  the  usual  way,  i.e. 

z  -  r  e“.  (2-5) 

where  r  is  the  modulus  and  the  argument  0  is  defined  as 

•  -  coe  -1  (u^Wl-32  ).  (2-6) 
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this  "pathological"  (computer  Jargon)  problem  is  to  be  desired  if  it  appears 
that  such  problems  might  arise  in  actual  practice. 


B.  Likelihood  of  Occurence. 

Before  attempting  to  attack  the  problem  of  clustered  roots*  it  is 
desirable  to  assess  the  possibility  of  their  actually  oc curing  in  practice. 
Digital  design  experience  shows  that  the  probability  of  system  roots  occur ing 
at  the  t  -  -1  singularity  is  unlikely*  unless  deliberately  placed  there  by  a 
digital  filter.  Since  the  system  designer  has  the  option  of  removing  such 
roots  by  altering  the  characteristics  of  his  digital  filter,  this  problem  can 
be  solved  readily.  However*  the  problem  is  not  so  readily  solved  if  the  roots 
are  clustered  at  the  z  •  +1  singularity. 


Suppose  that  a  system  plant  has  a  second  order  characteristic,  such 
as  one  would  find  typically  with  a  bending  mode  or  in  the  denominator  of  a 
transfer  function  (T.F.)  describing  the  characteristics  of  some  sensors  or 
actuators*  A  typical  transfer  function  might  be 

T.F.(s)  -  .  (2-9) 

(»fa)2  +00 2 

where  the  numerator  N(s)  may  have  any  farm  (in  this  discussion)*  An  example 
for  a  bending  mode  would  have  N(s)  represent  the  modal  gain,  nan  would  be  the 
modal  damping  ratio,  and  'V*>"  would  be  the  square  root  of  the  difference  of  the 
squares  of  the  modal  natural  frequency  (uj  n)  and  "a" ,  i.e* 

to  -  (to2  -  a2  )*  *  (2-10) 

n 

The  z-transfarm  of  the  transfer  function  of  Eq.  (2-9)  is  seen  to  be 

T.F.(z)  -  -2 - ,  (2-11) 

-  2ze  aTcos  *T  +  e  *aT 

where  T  represents  the  sampling  period  of  the  digital  system*  The  nature  of 
numerator  N(z)  depends  on  the  form  of  N(s).  The  important  thing  to  observe  is 
that  if  a  or  T  or  both,  (i.e*  if  their  product)  approach  zero,  then  the  denominator 
of  Eq.  (2-11)  can  approach  the  farm, 

D(z)  -  (z-1)2  .  (2-12) 

In  this  unfortunate  case,  the  roots  would  be  clustered  near  (or  on,  in  the 
limit)  the  singularity  at  z  -  +1 ,  giving  rise  to  the  problems  brought  forth 
in  Section  A  above.  Thus  it  unfortunately  appears  possible  for  the  clustering 
problem  to  actually  arise  in  a  physical  system.  The  probability  of  such  an 
occur ance  depends  on  the  actual  dynamic  nature  of  the  plant's  sensors,  actuators, 
and  dynamic  bending  properties.  Judicious  engineering  design  can  of  course 
obviate  such  problems  brought  forward  because  of  actuators  or  sensors,  if 
desired.  However,  changing  the  bending  characteristics  of  the  missile  plant 
can  be  difficult.  All  1s  not  lost,  however,  because  it  is  the  product  of 
a  and  T  in  Sq.  (2-11)  that  gives  rise  to  the  problem.  Hence,  again  the  designer 
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can  select  the  numerical  value  of  T  that  will  preclude  such  a  clustering 
problem.  If  other  design  constraints  conspire  against  the  digital  system 
designer,  he  might  desire  another  analytical  approach  that  might  solve  the 
clustering  problem.  Should  that  desire  materialize,  another  approach  is  proposed 
and  discussed  in  the  following  sections  of  this  report. 


) 

I 

I 

SECTION  m.  PROPOSED  SOLUTION. 

The  problem  investigated  in  Section  H  is  attacked  in  this  section 
by  applying  a  transformation  from  the  complex  z-dcmain  to  the  complex  w-dcmain. 

As  desired  in  Section  2.2  of  the  referenced  Technical  Requirement ,  it  is 
intended  that  this  technique  "mathematically  spread  the  system  poles  and  zero 
via  intermediate  complex  transform  techniques  to  eliminate  the  sensitivity." 

As  in  Section  H  (for  the  sake  of  simplicity) ,  one  may  again  assume  the 
dynamics  of  a  system  are  represented  by  the  three  roots  of  Eqs.  (2-1)  -  (2-2). 

The  usual  definition  used  in  mapping  from  the  z-dcmain  to  the  w-domain  is: 
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(1  +  *) 
(1  “  w) 


OD 


The  complex  variable  w  is  composed  of  real  and  imaginary  parts  which  may  be 
denoted  in  the  usual  manner. 


w  -  <TW  +  icow. 
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Applying  Eq.  (3-1)  to  root  z^  (2-1),  one  may  obtain  its  equivalent  location 
in  the  w-domain: 


(1  +  w) 

*  “  (1  -  w‘,  ‘  Z1  -  d  * 

If  one  solves  Eq.  (3-1)  far  w,  one  obtains  the  root  location 
(1  -  d) 

W.  as  •  • 

1  (1  +  d) 


03) 


(3-4) 


It  may  be  observed  in  passing  that  the  real  root  location  in  the  z-dcmain  has 
been  transformed  into  a  real  root  location  in  the  w-dcmain.  In  a  similar  manner 
the  equivalent  w-domain  expression  for  the  complex  conjugate  root  of  Eq.  (2-2) 
may  be  obtained: 


(1  +  w)  -  id 

8  "  (TTT)  "  z2 x 


+  iy 
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Sq.  (>-6).  The  three 
as  shown  in  Fig.  2. 


W. 


If  it  is  assumed  that  w^  #  w2  •  and  comprise  the  roots  of  the  same 
system  analyzed  in  Section  H,  then  the  characteristic  equation  may  be  written 
as 


C.E.(w)  -  (w-w^)  (w-w2)  (w-w^)  -  0.  (3“7) 

Alternatively,  the  characteristic  equation  may  be  determined  by  subsituting  the 
value  for  z  of  Eq.  (3-1)  into  C.E.(z)  of  Eq.  (2-?)  to  obtain 

C.E.(w)  -  (l+dt2x*--p2+2dxKi 

+  (3+d+2x-  _p2-2dx-3d  ^2)w2  +  (3-d-2x-J>2-2dxOd/2)w 

+  (l-d-2*f /2+2dx-d j>2)  -  0.  (3-8) 

One  observes  in  this  case  that ,  as  the  locations  of  the  three  roots  coalesce , 
i.e.  as  approaches  d,  one  obtains  (in  the  limit)  from  Eq.  (3-8)  the  expression 


lim  C.E.(w)  -lim  C.E.(w)  -  (l+d^w3  +  3(l+d)2(l-d)w2 

f-*d  J>*  d 

x->  d  +  3(l+d)(l-d)2w  +  (1-d)3  -  0. 

y— »  0 


09) 


If  d  approaches  unity  (as  before),  then 
lim  C.E. (w)  -  (2w)^  -  0, 
d-*l 


(3-10) 


and  the  three  roots  approach  (as  a  cluster)  the  origin  of  the  w-plane.  Now, 
however,  if  the  scale  of  the  system,  i.e.  —  in  the  case  of  Fig.  2--is  increased, 
success  is  met  in  spreading  apart  the  roots  because  the  vector  to  each  originates 
at  the  nearby  origin.  Recall,  in  the  z -domain,  increasing  the  scale  did  not 
substantially  spread  the  roots  apart  because  the  vectors  to  the  roots  went 
(approximately)  from  the  origin  to  the  z-+l  point. 


9 


amtpMMu:  -.^iMrcf^nKVRMi 


SECTION  IV.  COMPARISON. 

The  scaling  problem  and  proposed  solution  alluded  to  above  in  the  w-domain 
may  not  be  immediately  apparent.  Therefore  several  numerical  comparisons  are 
ma/fe  in  this  section.  Ten  runs  were  made  for  selected  values  of  roots  in  the 
z -domain.  The  corresponding  roots  in  the  w-dcraain  were  calculated  for  cote* 
parison.  The  results  were  obtained  using  the  program  of  the  Appendix  on  the 
Hewlett-Packard  9100B  Calculator.  They  are  summarized  in  the  following  Table. 

TABLE.  ROOT  LOCATIONS  IN  Z-  &  W-  DOMAINS 


Run 

i ^  »d 

z2“x*-iy 

W1 

W2 

mhi 

1 

0.9 

0.&H0.2 

-0.05 

-0. 104-10.12 

2 

0.96 

0.95+10.02 

-0.02 

-0.01+i0.03 

3 

0.98 

0.97+i0.02 

-0.010 

-0.015*i0.010 

k 

O.99 

0.9&*i0.01 

-0.005 

-0.010+10.005 

5 

0.98 

0.99*i0.01 

-0.010 

-0. 0054-10.005 

6 

0.996 

0.95*0.001 

-0.002 

-0.00251+15x10"^ 

7 

0.995 

0.996*10.001 

-0.00251 

-0.002 OO  iSXlO"4 

8 

0.999 

0.999*10.001 

-5X10"4 

-sxio^+isxio"4 

9 

0.9999 

0.999^10.0001 

-5X10"5 

-5I10'5+i5X10"5 

10 

0.99999 

0.99999*0.00001 

-5I10"6 

-5X1 0"^+ i5Xl 0”^ 

It  is  seen  that  the  distance  between  the  z  roots  and  the  singular  point  at 
z-+l  is  about  twice  the  distance  between  the  w  roots  and  the  origin  of  the 


It  is  readily  observed  that  increases  in  scaling  in  the  z-danain  do  not 
have  a  significant  effect  on  roots  clustered  about  the  z-+l  singularity. 
This  is  because  the  vectors  to  these  roots  eminate  from  origin.  The  same 
is  not  true  in  the  w-doroain.  Here  it  is  seen  that  large  increases  in  scale 
have  profound  effect  on  spreading  out  the  roots  clustered  about  the  origin. 
This  is  because  (in  this  case)  the  vectors  to  these  roots  eminate  from  the 
origin  to  the  (in  this  case)  nearby  root  locations. 
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SECTION  V.  CONCLUSION 


An  Investigation  has  been  made  to  investigate  the  potential  existence 
o if  a  clustering  of  system  poles  and  zeros  in  the  z-domain.  It  was  shorn  by 
example  in  Section  HA  how  such  a  problem  might  arise.  In  Section  113  the 
likelihood  of  such  an  occurance  was  assessed.  It  was  shown  that  such  a 
situation  could  occur,  but  that  there  are  design  techniques  to  prevent  it. 

If  these  fail,  an  alternate  approach  is  described  in  Section  H  whereby  the 
system  poles  and  zeros  are  transf armed  from  the  complex  z -dona in  to  the  complex 
w-dcmain.  It  then  is  possible  to  increase  the  scaling  (in  the  w-dcmain) 
with  significant effect  in  separating  the  clustered  poles  and  zeros,  thereby 
obviating  computational  difficulties  that  might  have  arisen  because  of  the 
clustering  in  the  z-domain. 
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